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Rigidity Percolation with g degrees of freedom per site is analyzed on randomly diluted Erdos- 
Renyi graphs with average connectivity 7, in the presence of a field h. In the (7, h) plane, the rigid 
and flexible phases are separated by a line of first-order transitions whose location is determined 
exactly. This line ends at a critical point with classical critical exponents. Analytic expressions are 
given for the densities riF of uncanceled degrees of freedom and 7r of redundant bonds. Upon crossing 
the coexistence line, jr and np are continuous, although their first derivatives are discontinuous. 
We extend, for the case of nonzero field, a recently proposed hypothesis, namely that the density of 
uncanceled degrees of freedom is a "free energy" for Rigidity Percolation. Analytic expressions are 
obtained for the energy, entropy, and specific heat. Some analogies with a liquid-vapor transition 
are discussed. 

Particularizing to zero field, we find that the existence of a. {g + l)-core is a necessary condition 
for rigidity percolation with g degrees of freedom. At the transition point 7c, Maxwell counting 
of degrees of freedom is exact on the rigid cluster and on the {g + l)-rigid-core, i.e. the average 
coordination of these subgraphs is exactly 2g, although 7c, the average coordination of the whole 
system, is smaller than 2g. 7c is found to converge to 2g for large g, i.e. in this limit Maxwell 
counting is exact globally as well. 

PACS numbers: 05.70.Fh, 64.60.Ak, 02.10.Ox 



I. INTRODUCTION 

Scalar Percolation(SP) P, 0] is a paradigm for the ge- 
ometric phase transition that takes place on an initially 
disconnected lattice of point-like sites, when the density p 
of present bonds is continuously increased. At the perco- 
lation point pc the system becomes connected on a macro- 
scopic scale. This may mean that transport can happen 
across the system (conductivity, fluid flow), or that the 
system becomes correlated on macroscopic scales, i.e. or- 
dered. It is in fact possible to describe magnetic transi- 
tions in terms of the percolation of properly defined clus- 
ters ^ 1^ . Because of the generality and simplicity of 
the concepts involved, this paradigm has found multiple 
applications in science 0, 0, IE H- most cases, the 
physical variables attached to sites are scalars, i.e. each 
site has one associated degree of freedom. 
A generalization of this paradigm considers the case in 
which there is more than one degree of freedom per site, 
and has been termed Rigidity Percolation(RP) j^, 0, 0, 
HI m m El m. in ^-RP, each site of a lattice has g 
degrees of freedom, and each present bond eliminates one 
relative degree of freedom. The most commonly invoked 
application of RP deals with the statics of structures. 
Consider for example the problem of bracing a frame- 
work in three dimensions, i.e. rigidly connecting a set of 
point-like joints by means of rotatable bars. Each joint 
has three translational degrees of freedom, g = 3, and 
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each rotatable bar fixes the distance between two nodes, 
thus providing one relative constraint. The question of 
whether a given set of bars is enough to rigidize a given 
structure constitutes a classical prob lem in applied math- 
ematics, that of Graph Rigidity [ia>li3jll3- statistical 
physics, bonds (bars) are randomly present with proba- 
bility p or absent with probability I ~p, and one asks for 
the typical rigid properties of such structures. Upon in- 
creasing the density p of present bonds, the system goes 
into the rigid phase, characterized by the existence of 
an extensive rigidly connected cluster. The RP problem 
becomes fully equivalent to SP when g — I. SP has a 
continuous transition in all dimensions 1] . In two dimen- 
sions the RP transition is continuous but in a different 
universality class 60] than SP. In the Mean Field(MF) 
limit, RP {g > 1) has a first-order transition [l^l20|. 
The situation is reminiscent of the Potts model, whose 
MF transition is continuous for q — 2 and discontinu- 
ous for g > 2 Potts models in the presence of 
a field, and their relation with percolation models, have 
been studied recently because of pos sible links with the 
deconfining transition in QCD 23,^^. In the presence of 
a field, and for large enough q, the Potts model has a line 
of first-order transitions ending at a critical point p5ll2^ . 
This critical point appears to always be in the Ising uni- 
versality class. 

Nonzero field values have been considered in scalar per- 
colation studies previously HHH IH SIM IH 111. 

A field may be introduced in percolation by allowing for 
the existence of "ghost bonds" which are present with 
probability h and connect sites directly to a solid back- 
ground (or to "infinity"). However, for any nonzero field 
there is no SP transition. RP is somewhat more interest- 
ing, as we will find out. In this work, RP is studied on 
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diluted random graphs of the Erdos-Renyi type, with av- 
erage connectivity 7, in the presence of a field h. Unlike 
SP, mean- field RP has, in the presence of a field, a line 
of first-order phase transitions. This line ends at a crit- 
ical point with classical critical indices: a = 0, = 1/2, 
7 = 1 and 5 — 

Some of the analysis in zero field is relevant for the related 
problem of Bootstrap Percolation (BP) 34, 35, 36, 37], 
also known as A;-core 38] in the field of Graph The- 
ory (3i,|4l|43. In fc-BP j3Jj, all sites with less than k 
neighbors are iteratiyely culled. What remains, if some- 
thing, is the fc-core js^; a subgraph were all sites have 
k or more neighbors. In g-RP a site needs at least g 
bonds in order to be attached to a rigid cluster. Thus 
the "infinite" rigid cluster is a subset of the g-core. In 
Section III Bl we will see that an even stronger condi- 
tion exists for rigidity: g-RP requires the existence of 
d. [g + l)-core 61]. 

Our approach starts by deriving an equation of state for 
the "order parameter" R{j,h), the probability that a 
randomly chosen site belongs to the rigid cluster, as a 
function of 7, the average number of bonds impinging on 
a site, and h = jH, where H is the applied field. We will 
call H or h indistinctly the "field" variable. The equation 
of state, as is customary in these cases, is found to ac- 
cept multiple solutions. Stability analysis is not enough 
to single out a unique solution. Of central importance 
in order to lift this multiplicity are 7^(7, /i), the average 
number of "redundant" bonds per site (see later), and 
np, the average number of uncanceled degrees of free- 
dom per site. Their relevance resides in the fact that 
they must be continuous functions of 7. Requiring that 
np (or, equivalently, 7^) be continuous is enough to iden- 
tify the physically correct solution. 

This work is similar in spirit to pre vious treatments of 
RP in zero field on Bethe lattices [23 , i-e networks where 
each site has exactly z randomly chosen neighbors. We 
consider the effect of an external field, and particularize 
to random graphs of the Erdos-Renyi type |42l ] . The in- 
troduction of a field appears to be much more tractable 
analytically on Erdos-Renyi graphs than on Bethe lat- 
tices, and this is the main reason why most of the results 
presented here are for Erdos-Renyi graphs. For these, 
we are able to derive analytic expression for the densities 
rii;- of uncanceled degrees of freedom and 7^ of redun- 
dant bonds for arbitrary values of the field. However in 
zero field calculations are also straightforward for Bethe 
lattices, and allow us to obtain, in an entirely analytic 
fashion, some of the results obtained by numerical inte- 
gration in Ref . 20] . 

It is known that SP and Potts models are particu- 
lar cases of a more general Fortuin-Kasteleyn Random- 
Cluster(RC) model, defined by a continuous parameter 
q 0! Q]' SP can be obtained as the g — > 1 hmit of 
this Random Cluster model and, in this limit, the loga- 
rithm of the partition function coincides with the average 
number of connected clusters. It seems possible that a 
similar mapping might exist for RP as well, although it 



has not been found up to now. However it has been pro- 
posed ,2^ ^] that the number of uncanceled degrees of 
freedom np is a good "free energy" candidate for RP in 
zero field. When g = 1, each connected cluster has one 
uncanceled degree of freedom, so both definitions coin- 
cide in this limit. We explore on this idea further in this 
work. 

It is possible to establish a pedagogical parallel between 
the RP transition and a Condensation transition. One 
identifies the coordination parameter 7 with an inverse 
temperature /3; the negative of the order parameter R 
plays the role of the fluid volume V (it is also possible 
to identify R with p, the fluid density) [g^] , and the field 
H = h/j is the fluid pressure P. Within this analogy, 
it results natural to argue that requiring the continuity 
of np in order to identify the physically correct solution 
R in Rigidity Percolation is equivalent to requiring the 
continuity of the free energy (giving rise to the Maxwell 
construction) in the Statistical Mechanics treatment of 
condensation transitions at the MF level [46^ . 
It is shown in this work that the idea of identifying np 
with a free energy for RP leads to consistent results in 
the presence of a field as well. In Section IVTI it is shown 
that: 1) the order parameter R is obtained as a deriva- 
tive of the free energy with respect to the field, 2) the 
condition of stability (Section IIII C|) may be related to 
the positivity of a suitable second derivative of the free 
energy, and, 3) the continuity of np can be cast exactly 
in the form of Maxwell's Rule of Equal Areas on the RP 
equivalent of a P-V diagram (which is the H-R diagram, 
see Fig. [7||. 

This work is organized as follows: The system under con- 
sideration is defined in Section III Al and its equation of 
state is derived in Section IIIBI establishing its connec- 
tions with Bootstrap Percolation. In Section lll CI the field 
is introduced. Section ITTll starts the analysis of solutions 
of the general equation of state, discussing stability and 
the existence of a critical point. The concept of redun- 
dant constraints is introduced in Section llV Al and their 
density is calculated in Section llVBI This result is used 
in Section to determine the value jc where the first- 
order transition takes place. In Section lVBI the counting 
of constraints is done on the (g+l)-Tigid core in zero field, 
both for Erdos-Renyi graphs and for Bethe lattices. Sec- 
tion discusses several consequences of identifying the 
density of uncanceled degrees of freedom with a free en- 
ergy, and Section IVlTl cont ains a discussion of the results. 



II. SETUP 
A. Randomly diluted graphs 

We consider graphs made of sites (or "nodes" ) where 
each of the N{N—l)/2 pairs of distinct sites is connected 
by a bond (or "edge" ) independently with probability p. 
This defines ^3] an Erdos-Renyi graph with average co- 
ordination number 7 = p{N — 1). In this work 7 will be 
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taken to be of order one. For large N, a site of this graph 
is connected to k other sites with poissonian probabiHty 
Pfe(7) = e-^7Vfc!. 

As appropriate for Rigidity Percolation, each node of this 
graph is regarded as a "body" with g degrees of free- 
dom. For example, rigid bodies in d dimensions have 
d translational degrees of freedom plus d{d — l)/2 rota- 
tional degrees of freedom, for a total of g = d{d + l)/2. 
In this work, g is taken to be an arbitrary integer. Each 
present bond represents one constraint, that is, removes 
one degree of freedom. In order to have a physical repre- 
sentation in mind, a bond can be though of as a rotatable 
bar that fixes the distance between two arbitrary points 
belonging to the bodies it connects 0,E3|- For the class 
of graphs that we consider, at most one bond is allowed 
for each pair of nodes. 

As the number of present bonds grows, parts of this graph 
will become "rigidly connected" . Rigid connectivity of a 
subgraph means (this is our definition) that the total 
number of degrees of freedom in the system cannot be 
further reduced by additional bonds connected between 
nodes of this subgraph. Such subgraphs are customarily 
called "rigid clusters" . A rigid cluster has no internal de- 
grees of freedom left. When the coordination 7 is large 
enough, the largest rigid cluster encloses a finite fraction 
of the system's sites, and rigidity is said to percolate. For 
low 7 the system is in the flexible, or floppy, phase, and 
there are no extensive rigid clusters. 
A rough estimate of the threshold for the appearance 
of an extensive rigid cluster can be obtained by equat- 
ing the average number of constraints per site, which 
is 7/2, to g. This results in the so called Maxwell es- 
timate m El El, 7^<»^'"«" = 2g. This estimate 
becomes exact in the g >> 1 limit fSection IV Bp . 



B. The Equation of State in zero field 

Let us now define C to be the largest rigid cluster in the 
graph, if one exists. A site in C will be said to be a "rigid 
site". Let R be the probability that a randomly chosen 
node i be in C In order for i to be in C it is necessary and 
sufficient that it be connected to g or more other nodes 
j in C. Our definition of rigidity is thus recursive at this 
stage. 

A site j is rigid and connected to i with probability pR. 
Therefore i has exactly k rigid neighbors with probabil- 
ity (\^) ipR)''{l - pR)^~^-^. For N large and defining 
X — jR, the probability Pk to have exactly k rigid neigh- 
bors may be written as 

Pk{x) = e-^x'^/kl, (1) 

showing that the number k of rigid neighbors of a ran- 
domly chosen site is a poissonian variable with average 
X = 7i? = X]fc=o ^-^k (x) ■ Since a site must have g or more 
rigid neighbors in order to be itself rigid, we conclude that 



R satisfies the self-consistent equation: 

00 

R = ^Pfe(7^) = Gg(7i?), (2) 

fc=g 

where we have defined 

00 

GM = J2 Pk{.x). (3) 

Bootstrap Percolation, or k-core 

Let us now brie fly d iscuss the related problem of Boot- 
strap Percolation Is^l or fc-core[38|. We want to assess 
the probability BP{g + 1) that a randomly chosen site i 
be part of the (5 -I- l)-core. Assume this is the case. By 
following one of its links, a neighbor j is reached, j must 
have at least g other neighbors. We call these neigh- 
bors of j, other than the site i from which we arrived 
at it, the "outgoing" neighbors of j. Each of these in 
turn must have g or more outgoing neighbors, and so on. 
More formally, let us define the property of y-outgoing- 
connectedness (g-OC) in the following (recursive) way: 
A site is g-OC if g or more of its outgoing neighbors also 
are. 

On graphs of the Erdos-Renyi type, the probability to 
have k outgoing neighbors is the same as that to have k 
neighbors altogether, i.e. Pk{l) (Eq. (Ql). This is so be- 
cause links are independently present. Letting R be the 
probability that a site be g-OC, by the same reasoning 
as in the previous Section we conclude that a random 
site is connected to exactly k outgoing g-OG neighbors 
with probability Pk{'jR)- Thus R satisfies the same self- 
consistent equation ^ as R does in the RP problem with 
g degrees of freedom, i.e. 

R=Gg{jR). (4) 

Notice that, by requiring that each site in a tree have g 
or more g-OC neighbors we ensure that all sites, except 
perhaps for the top one, have {g + 1) or more neighbors. 
The probability BP{g -f 1) that the top site itself has 
g + 1 or more neighbors (which are g-OC) is then given 

by 

BP{g + l) = Gg+^{-iR) = R-Pg{-iR). (5) 

This expression gives the density of the {g + 1)-BP infi- 
nite clus ter (or (q + l)-core) at the point where it first 
appears [3J, |33, EJ . So one must first solve Q in order 
to obtain ii as a function of 7, and then use jSJ to find 
BP{g + 1). Numerical results |63| show that is exact 
for large N . 

We see that the [g -\- l)-core density PB{g + 1) is some- 
what smaller than R whenever i? > 0, while R in turn 
satisfies the same equation as the density R of rigid sites 
in .g-RP. Later in Section IV Bl we will see that in zero 
field, whenever there is a g-rigid cluster, it contains as a 
subset the {g -\- l)-core. 
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C. Equation of State in the presence of a ghost 
field 

We now introduce a "ghost field" H that couples to 
the order parameter R, in the following way |32j| : in 
addition to the "normal" bonds of our graph, a num- 
ber N^H — Nh of "ghost-bonds" are assumed to ex- 
ist, each connecting one site to a unique rigid back- 
ground |64| . Each ghost bond provides one constraint, 
i.e. removes one degree of freedom. Multiple occupa- 
tion is allowed, so that a random site is connected to 
the background by n ghost bonds with poissonian prob- 
ability Pn{h) = e-''h"/n\. 

It is easy to see that all extensive rigid clusters are rigidly 
connected to the background. From now on, a site is said 
to be rigid if it is rigidly connected to the background, ei- 
ther directly through ghost bonds, or indirectly through 
rigid neighbors. 

If a site has n > g ghost edges connecting it to the back- 
ground, then it is rigid. Otherwise if n < g, it is rigid 
if it has, in addition to these n ghost bonds, {g — n) or 
more rigid neighbors. Thus using ^ we may write, in 
the presence of a field h — H'y, 

oo 5~ 1 oo 

n—g 71—0 k—g — n 

After a simple resummation, this expression reads 

R^Gg{^R + h), (7) 

with Gg given by © . Eq. (jZj) generalizes ^ in the pres- 
ence of a field, and is the equation of state for our prob- 
lem. 

The simplicity of Q is one of the reasons leading us to 
study this particular (Poissonian) definition of the field. 
Other field definitions, like for example assuming that 
each site is rigidly attached to the background with prob- 
abilit y h , or other random graph structures like Bethe lat- 
tices |23| are also tractable with the methods used here, 
but the algebra becomes more complicated. 
Clearly y — ^{H -I- R) plays the role of a "Weiss field" 
in the MF equation for a ferromagnet. By analogy we 
may thus identify 7 as the inverse temperature, and H 
as the magnetic field. It is illustrative to take notice of 
the similarities between ((TJ and other MF equations. For 
the Potts model 

m = eqil3m + h), (8) 

where (3 is the inverse temperature, m is the mag- 
netization, h = (3H is the external field, and 

egiy) = {ey-i)/iey + {q-i)). 

Notice that when g — I (Scalar Percolation), ||7J) gives 
Gi{y) = 1 — , the same as ^ ior q = 1. This is of 
course just a consequence of the known equivalence be- 
tween Scalar Percolation and the q ^ 1 limit of the Potts 
modellillil. 



However, the archetypal example of a first-order transi- 
tion with a two-parameter phase space is the condensa- 
tion transition , described at the Mean-Field level by 
the Van der Waals equation. The reduced form of the 
Van der Waals equation can be written as 

\^V(fi{\ + P)), (9) 

where r(y) = ily"^ /{y -|- 8)^, P is the pressure, /3 is the 
inverse temperature and A = with p the number den- 
sity, r has the same general features as G and O, namely 
it starts at zero, grows sharply, and saturates for large 
values of its argument. Therefore all three systems give 
rise to the same phenomenology, including of course shar- 
ing the same (classical) critical indices at their critical 
points. 

Within an analogy with a condensation transition, in the 
RP problem 7 plays the role of an inverse temperature, 
while H — h/"f plays the role of a pressure. 



III. ANALYSIS OF SOLUTIONS 

In this Section, the solutions of our equation of state 
(|7|l are discussed. In order to obtain i?(7, h) numerically 
for given 7 and h, one might for example iterate Q until 
a desired numerical accuracy is reached. This procedure 
was used in Refs 0, |23| . However in this work the fol- 
lowing alternative procedure was preferred: given h > 
fixed, and for a sequence of values of ?/ > h, evaluate 

R = Ggiy) 

(10) 

7 = {y-h)/GM 

thus defining i?(7, h) implicitly. In this way one obtains 
the results displayed in Fig. ^ This procedure allows us 
to obtain all solutions of iQ, while the iterative proce- 
dure mentioned above only provides the stable branches. 

In zero field, i? = is a solution V7 and Vg. If 7 is large 
enough, nontrivial solutions R> exist as well. If /i > 
one has on the other hand that i? > V7. 
In the SP case {g = 1), Gi{x) ~ 1 — exp{—x). Thus if 
h — Si non-zero solution for R first appears at 7c = 1. 
Above 7c one has i? ~ (7 — 1). In this case there is a 
continuous transition at 7 = 1, a well known result for 
scalar percolation on random graphs "i^]. If ft, > there 
is no transition (Fig.^l. 

A richer behavior is found for the RP case {g > 1), as 
shown in Fig. ^ for g = 2. In this case 7, as given by 
Eq. may no longer be a monotonous function of y, 
and, for this reason, R becomes a multivalued function of 
7. This allows for the existence of a first-order transition. 
The physical considerations which lead to the identifica- 
tion of the correct transition point will be discussed in 
Sections Imn and IVl 
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FIG. 1: The density R of "rigid sites" as given by Eq. Q, for 
g = 1 (scalar percolation, top) and g = 2 (rigidity percolation, 
bottom). The field takes the values (from right to left): h — 
0.0, 0.1, 0.20, . . . , 1.0. For g = 2 the critical field is he = 0.287. 
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FIG. 2: For all g > 1, the order parameter R as given by 
l|10|l becomes multivalued for small values of the field h. The 
branch joining the spinodals 7± (dashed) has dR/dy < and 
is thus unstable. The uppermost and lowermost branches 
(solid) are stable. In this example, g = 2 and h — 0.15. For 
zero field the spinodal 7_ goes to infinity, and the branch 
< 7 < 7- collapses onto the R — solution, which is stable 
for all 7. 



B. Critical Point 

When h takes a critical value h*, the spinodals x^. 
coalesce onto an inflexion point where dGg{x + h)/dx 
equals Gg/x (since Eq. (|ll|l holds). The critical point 
{/i*, 7*, i?*} is thus defined by 



A. Spinodal Points 

The condition that 7 as given by Eq. pU|l be stationary 
in X reads 



G„ 



d_Gg 
dx 



= 0, 



(11) 



and has two solutions x'!^{g,h) for all h smaller than a 
critical value h*{g). These two solutions in turn define 
7|., which are turning points for R{j, h) (See Fig.|2). In 
the interval 7^ < 7 < 7^, R{j,h) has three solutions, 
two of which are stable as we show next. Therefore a 
discontinuous 7-driven transition takes place when h < 
h*(g). Right at h — h* the transition becomes continuous 

"(g), i? is a smooth 



(See Section IIITB(I . while for h > h* 
function of 7 and no transition occurs. 

The thresholds 7^. can be identified as spinodal-points, 
as first discussed in ■ The true rigidity percolation 
transition happens at a value 7c (5, h) that lies in between 
the spinodals, and which we analytically determine later 
m Section El 



dx'^ 
dG„ 



dx 





Gg{x*+h*) _ J_ 

X* 7* 



(12a) 
(12b) 



For h > h* there is no phase transition. Using Eq. 0, 
and defining y = x + h, these two equations can be solved 
exactly. The critical point turns out to be 



r 2/* = .9-1 

7* = {P3-i(.g-l)}-i=eff-i(g-l)!/(g-l)f-i 

R* = G,(5-l) = ei-sEr=,(5-l)Vfc! 

L h* = y* -X* =g-l- 7*i?* 

(13) 

For the scalar case (g = 1) this gives h* = R* = and 
7* = 1 in accordance with the fact j43] that for zero field 
there is a continuous transition at 7c = 1. 
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For g = 2 one finds 

' 7* = e 

< R* = (e-2)/e (14) 

^ h* = 3-e. 

For large g, and approximating n! « (n/e)"(27rn)^/^ one 
sees that 7* cx g^/^, while h* ~ g and thus H* ^ g^l"^. 
This means that in the limit 17 >> 1 most constraints are 
field-constraints at the critical point. 



1. Scalar Percolation (g = 1) 

In zero field, i? = is a solution of Eq. {Tj), thus a 
fixed point of Eq. IIH), V5. If 5 = 1, stability (Eq. (|T7|)) 
requires that 7P0 = 'je'"'^' < 1. Thus the trivial solution 
R = becomes unstable for 7 > 1, where the nontrivial 
solution (stable because dR/dj > 0. See Fig. ^ first ap- 
pears. This situation is typical of continuous transitions; 
the ordered solution appears exactly at the point where 
the paramagnetic solution becomes unstable. 
When h > there is only one solution for Eq. {Tj), and it 
is stable for aU 7 since dR/dj > (Eq. l(TS|) '). 



1. Critical Indices 

The RP transition on random graphs is similar to other 
MF transitions with first-order lines, as discussed in Sec- 
tion Thus it can be concluded that RP must have 
classical critical indices: (3 — 1/2, S = 3 and 7 = 1. 
For completeness we show that this is indeed the case by 
deriving the critical indices briefiy in Appendix IbI 



C. Stability analysis 

Stability can be analyzed if Eq. |(7J is interpreted as a 
recursion relation, 



R 



t+i 



G{^Rt + h), t = 0,l,. 



(15) 



for fixed 7 and h. Assume R{^, h) is a fixed point of 
Eq. (fT^ . i.e. R — G{^R, h). This fixed point is stable if 



dG(nR + h) 



OR 







7,/i 





dG{x + h) 



dx 



x—'yR 



< 1. 



(16) 



Since G'{y) — Pg^x{y) we conclude that the stability 
condition reads 



iPg-liv) < 1, 



(17) 



with y = X + h as defined previously. 

A more useful form results if one notices that, if i? > 0, 

Eq. H16|) is equivalent to requiring that dR/dj > 0. In 

fact, 



dR 

dj 



dG dx ^1 d 

dx c?7 c?7 



{iR) = G' [r + -i 



dR 

dj 



^(1-7G') =G'i? 
07 



(18) 



and since G" > Vx > 0, we conclude that a nontrivial 
fixed point i?(7, h) > oi Eq. H15|) is stable if and only if 
dR/d-y > 0. 

This condition has a simple physical meaning: increas- 
ing the average connectivity 7 should not decrease the 
rigid density R. A similar stability condition holds for 
fiuids, namely that the coefficient of thermal expansion 
be positive. 



2. Rigidity Percolation (g > I) 

For g > I the situation is more interesting, as several 
stable solutions of Eq. iQ can coexist. U h < h* , Eq. O 
has three solutions (see Fig.[2Jl in the range 7! > 7 > 7iJ_. 
The branch joining the two spinodal points can be dis- 
carded since it is unstable (because dR/dj < 0). However 
the other two branches (solid lines in Fig. O are stable 
so that there is coexistence of two stable solutions for 
7* > 7 > 7^. Thus the system undergoes a 7-driven 
first- order transition somewhere between the spinodals. 
The precise point 7c at which the system switches from 
one stable solution to the other is uniquely defined by 
a continuity requirement, as we elaborate later in Sec- 
tion El 

The h = case is similar, however in this case the lowest 
stable branch collapses onto the trivial solution R — 0, 
while the spinodal 7_ goes to infinity. 



IV. ZERO MODES, REDUNDANT AND 
OVERCONSTRAINED BONDS 

In this Section the counting of uncanceled degrees of 
freedom and the useful notions of redundant and over- 
constrained bonds are discussed. These will be of central 
importance in our subsequent treatment of the RP prob- 
lem. 



A. Definitions 

We consider a graph made of N sites, each with g de- 
grees of freedom, and for the moment assume that h = 0, 
i.e. there are no ghost-field constraints. Each present 
bond removes one degree of freedom from the system, 
unless it is a redundant bond. A redundant bond (or 
constraint) is one that links two nodes which were al- 
ready rigidly connected, as for example nodes i and j in 
Fig. El The addition of a redundant constraint does not 
reduce the number of degrees of freedom in the system. 
Thus the balance of degrees of freedom reads 



Nf = Ng - E + Br, 



(19) 
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where A'"^' is the number of uncanceled degrees of free- 
dom [65|. E is the number of bonds in the graph, and Br 
is the number of redundant bonds. 

For an Erdos-Renyi graph one has (E) — 'yN/2. Defining 
{Br) — jrN/2, the average density np{'^) of zero modes 
per site is written 

nF{l)^^-^=g-\b-lr{l)]. (20) 
Let us assume that 5 is a redundant bond. This means 




overconstrained bonds, on the other hand, needs not be 
continuous. 

Let us now discuss how the balance of degrees of freedom 
is modified by ghost-field bonds f Section III C|) . If a site 
has n < g present ghost edges connecting it to the rigid 
background, then it contributes with (g — n) degrees of 
freedom to Eq. H19|) . If on the other hand it is connected 
to the background hy n > g ghost bonds, then this site 
has no degrees of freedom left to contribute to Eq. ifT^ . 
The effective number ges{h) of remaining degrees of free- 
dom contributed by an average site is thus 

a 

gcs{h) = ^{g-n)Pn{h) 

= g{l-Gg+^{h))-h{l-Gg{h)), (21) 

where Gg is defined by (|3J). Taking this result into ac- 
count, in the presence of a field h we now have that 



FIG. 3: In this two dimensional example, each node (circle) 
represents a point that has two positional degrees of freedom, 
while each bond fixes the distance between two points and 
thus provides one constraint. The graph on the left still has 
three remaining degrees of freedom: two translations and one 
rotation. Bond ij (right) is a redundant bond, because it con- 
nects two nodes which were already rigidly connected. After 
adding bond ij, the bonds that provided the rigid connec- 
tion between i and j become overconstrained (thick lines). 
Any one of these can be removed without altering the num- 
ber of remaining degrees of freedom, which still equals three. 
The graph on the right has one redundancy, but 12 overcon- 
strained bonds. 

that its two end-nodes are rigidly connected even if b is 
removed. Let Bb be the subset of bonds that, in the 
absence of &, provide rigidity to its two end-nodes. Af- 
ter adding b, any of the bonds in {Bb + b} (thick lines in 
Fig. 01 can be removed without altering the total number 
of degrees of freedom Np. The subset of bonds {Bb + b} 
is said to be overconstrained. From a physical point of 
view, the overconstrained bonds may be defined as those 
that carry an internal stress because of the addition of a 
redundant bond b that has a length-mismatch. 
Please notice the important difference between the num- 
ber of redundancies and the number of overconstrained 
bonds: when adding a redundant bond (like ij in Fig.EJ 
the number of redundancies Br always increases by ex- 
actly one. However the number of overconstrained bonds 
may increase by more than one. In the example shown 
in Fig. 121 the number of overconstrained bonds increases 
by twelve. 

We then conclude that, when adding any bond to a 
graph, the number of redundancies Br will either in- 
crease by one (if the chosen sites were rigidly connected) 
or stay unchanged (if not), in which case the number of 
zero modes Np will either stay unchanged or decrease 
by one. This implies that the densities 7^ of redundant 
bonds and np oi zero modes must be continuous func- 
tions of the density 7 of present bonds. The density of 



np{-f,h) 



(22) 



Notice that np and 7^ are trivially related through the 
addition of a known function of 7 and ti. Therefore the 
knowledge of either np or 7^ provides exactly the same 
physical information about the system. 
The importance of H22|l is twofold: in the first place, as 
discussed above, np and 7^ must be continuous functions 
of 7. A similar reasoning allows one to conclude that they 
must also be continuous functions of h. The continuity 
of "fr and is a key property, and will be used later in 
Section El to select the physically correct solution of O 
whenever indeterminacies arise. Secondly, in Section I VII 
it will be argued that np can be identified with the loga- 
rithm of the partition function for the RP problem, and 
the consequences of such identification will be discussed. 
Although this identification is not necessary for solving 
the RP problem, it provides interesting additional insight 
into this problem insofar it helps making a link with ther- 
modynamics. 



B. Calculation of 7^ 

We now show how the density 7^(7,^1) of redundant 
bonds is calculated. Obviously when 7 = there are no 
redundant bonds, i.e. 7r(7 = 0, /i) = 0, so we can write 



djr_ 

dj 



(23) 



where the integral is done along a path of constant h. In 
Appendix^ we show [b^ that 



7- 



d{Br) 



97 



(Bov), 



(24) 



where {Bov)-f is the total number of overconstrained 
bonds in B. Consider two randomly chosen sites i and 
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j. As discussed in Section IIVI a bond 6.^ is overcon- 
strained if i and j are rigidly connected to the back- 
ground, even in the absence of this bond. On a ran- 
dom graph, this happens with probabihty R^. Thus 
a bond is present and overconstrained with probabihty 
{■-y/N)E?, and therefore the average number of overcon- 
strained bonds is {Bov) = We can now write 
Eq. 1(23 as 



so that now Eq. reads 



jr{l,h) = / R\^,h)dj. 



(25) 



(26) 



Notice that, while (|23|l and H24|l are exact for arbitrary 
graphs, H25|) only holds under the assumption that two 
neighboring sites i and j are independently rigid with 
probability R. Thus in deriving (|26|) . correlations have 
been ignored. This is correct on random graphs of the 
type considered here, but fails on finite dimensional lat- 
tices. 

It is not immediately obvious how (|26|l can be integrated 
since we do not have an explicit expression for i? as a 
function of 7. However using x — ^R — y — h one can 
write R^d'^ — Rdx — xdR = {R—{y — h)dR/dy)dy, where 
R = Gg{y) as given by Q. After integrating by parts. 



y—x-\-h 



7. - |2y Ggiy)dy ^ {y - h)Gg{y) 

Using Eq. H27|l can be integrated to give 

7. = {{y + h-2g)GM + Wy)}\lZT'' 
= {y + h)Gg{y)~2hGg{h) + 2g[Gg+i{h) 



(27) 



y=h 



Gn 



-i(y)] 

(28) 



This expression enables us to calculate the density of re- 
dundant bonds analytically. In previous work on Bethe 
lattices [l^, ()26|l was solved in zero field by numerical 
integration, using for this purpose the values of R that 
result from iteration of Although the algebra be- 
comes slightly more complicate, the approach used here 
to calculate 7,. is applicable to Bethe lattices as well. 
Later in Section FV Bl we present some results for Bethe 
lattices. 



1. The density of zero modes 

Using O, lEU, (US, and I^HIl we can write 

nF{-i, h) ^l{R^-l) + {h- g){R - 1) + gPg{y) (29) 

where y = x + h — ^R+ h = ^{R + H). Because of (|^ . 
continuity of 7^ implies that of n^?. 



2. A eonsistency check 

In the limit 7 ^ A^, all bonds are present and there are 
no remaining degrees of freedom, i.e. np ~^ 0. It is easy 
to verify that Eq. I|29|) indeed satisfies this condition. 
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FIG. 4: The density of redundant bonds 7r (thick dashed) 
describes a Maxwell loop (dot-dashed) whose crossing point 
defines the location of the true transition. After discarding 
the unphysical loop, the density R of rigid sites (thick solid) 
jumps at the transition between A and B, but 7^ is continuous. 
Shown is an example with g = 10 and h = 5. 



V. LOCATION OF THE TRANSITION POINT 

A. The general case 

The value 7c(/i) where the transition takes place is uni- 
vocally determined by the requirement that 7^ be a con- 
tinuous function of 7. In order to see how this works, 
imagine that for a given h < h* one lets h < y < 00 and 
calculates jr{y, h) (Eq. (|2Hl) and 7(1/, h) (Eq. COJ)), thus 
obtaining a plot of 7^. vs 7 (Fig. ^ . In doing so one finds 
that 7r describes a loop (sometimes called a "Maxwell 
loop" in Statistical Mechanics). 

The existence of a crossing point implies that two values 
yA and ys exist such that 



{{v + h-2g)Gg{y) + 2gPg{y)}\ 



y=VB 



0, 



(30) 



where yA.B — Jc{h)RA,B + h. Eq. cannot be solved 
analytically in general, however it is easily solved numer- 
ically (by iteration) for each value of h and g, to obtain 
lc{g,h). In this fashion the location of the coexistence 
line can be calculated to arbitrary precision (Fig.O. At 
the transition point "fc{h) the order parameter R jumps 
(Fig. 1^) , but the density of redundant bonds 7^ is con- 
tinuous (Fig.Eb). 
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FIG. 5: Coexistence lines (lines of first-order phase 
transitions-dashed) in the space defined by the two inten- 
sive parameters of Rigidity Percolation: the "field" H = h/'y 
and the "temperature" I/7, for several values of the number 
g of degrees of freedom per site. The critical point, defined 
by Cnji is indicated with an asterisk. For all g > 2, i.e. for 
Rigidity Percolation, this critical point has classical indices 
a = 0, P = 1/2, 7 = 1 and S — 3. The case g — 1 (Scalar 
Percolation, not shown) has its critical point at 7 = 1 and 
H — 0, and belongs to a different universality class, with 
a = -1, /3 = 1, 7 = 1 and 5 = 2. 



B. The case of zero field. Maxwell counting on the 
{g + l)-rigid-core 

For zero field one has that i? = for 7 < 7c. Eq. H26() 
then implies that 7^ = for 7 < 7c- At the transition 
point 7c, a rigid cluster suddenly appears that contains 
a fraction Rc of the graph. In zero field, Eq. H3U|I reads 



(xc - 2g)Rc + 2gPg{x,) = 0, 



(31) 




0.0 1.0 2.0 



3.0 

Y 



4.0 5.0 6.0 



where all quantities are evaluated in the rigid phase {R 
and y are zero in the floppy phase). 
Even this simpler equation does not admit analytic treat- 
ment. However, very precise solutions can be obtained 
numerically as described already. Table contains some 
examples. 

Please notice that 7c is smaller than the value 2g needed 
for global balance of constraints and degrees of freedom. 
However, all degrees of freedom belonging to sites on 
the rigid cluster (which appears suddenly at the transi- 
tion point) are by definition canceled out by constraints. 
Moreover, since for zero field the number of redundant 
constraints per site is exactly zero at the transition (See 
Fig. 0) , one concludes that the average number of bonds 
per site, for sites on the rigid cluster, must be exactly 
2g at 7c. The conclusion is then that sites on the rigid 
cluster have more bonds than average. 
This observation may be formalized in the following way. 
Using jR = x,R = J2Zg Pk{x) with Pk{x) = e-'^x^/k\, 
and xPk{x) — (fc + l)Pfc_|_i(x), condition (|31|l can be writ- 



FIG. 6: a) The final values for the density of rigid sites 
R (thick lines). The transition point 'ydh), where R un- 
dergoes a discontinuous jump, is determined by the physi- 
cal requirement that the density of redundant bonds 7^ (b)) 
be continuous. In this example g = 2 and (from right to 
left): h = 0.0, 0.1, 0.20, . . . , 1.0. For g = 2 the critical field is 
he = 0.287. 



ten 



^ {2g-k)Pkix)^0, 



(32) 



25- 



(33) 



which in turn implies 

T,'^=g+i Pkix) 

This sum counts the number of rigid neighbors for sites 
on the (g + l)-rigid-core, that is, the subgraph of the rigid 
cluster that has minimum coordination {g + 1). Eq. 133|) 
means that the {g + l)-rigid-core has an exact balance of 
degrees of freedom. 

Clearly, the rigid cluster may also contain sites not in 
the {g + l)-rigid-core. When these are considered, it is 
possible to see that the balance of constraints and degrees 
of freedom is still respected. 
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ic{g) 


discontinuity Rc 


2 


3.58804747296539 . . . 


0.74915378510250 . . . 


3 


5.7549256115462 . . . 


0.8812398370507 . . . 


4 


7.84295819428849 . . . 


0.933538491167168 . . . 


5 


9.89551361859910 . . . 


0.959638549916489 . . . 


6 


11.9288724790262765 . 


. 0.974275720961391 . . . 


—* oo 


2ff 


1.0 



TABLE L The Rigidity Percolation threshold coordination 7c, 
and the jump in the rigid cluster density Rc on Erdos-Renyi 
random graphs, obtained from solving 13111 for several values 
of g, the number of degrees of freedom per site. 



The case of large g 

An analytic solution of Eq. H31|l is possible (though 
somewhat trivial) in the g » 1 limit. In this case Pg — > 
and i?c — * Ij i-e. there is a jump from R = Q to R— 1 
at the transition point. Plugging these observations into 
(|S1Tl one gets, for zero field. 



7c 



2.9 



.9 >> 1- 



(34) 



This means that the rigid transition happens exactly at 
the point where global balance between constraints and 
degrees of freedom is attained. In other words, at 7c 
Maxwell counting is exact on the rigid cluster for all g, 
while for g >> 1 it is globally exact as well. 



Bethe lattices 

A graph where each site has exactly z randomly chosen 
neighbors is topologically equivalent to a Bethe lattice, 
except for finite-size corrections. When this system is 
randomly diluted, one obtains what is called "random 
bond model" in Ref. in which each of the z bonds 
on a site is present with probability p. We will only 
state the final results for this system. In zero field the 
transition point Pc is determined by 



E 



(35) 



where x=pT, Qf\x) = (: 
the recursive equation | T9| 



)x^(l 



3=a 



and T satisfies 



(36) 



The transition condition H35|l gives rise to a polynomial 
equation which is easily solved for small values of g and 
z. For g = 2 and z = 5 the corresponding quadratic 
equation results in pc = 0.83484234. For g = 2 and 
z = 6 a cubic equation is obtained, and its solution is 



Pc = 0.656511134. This last value is consistent with but 
more precise than pc = 0.656 as obtained in previous 
work nji u sinEf numerically exact matching algorithms 
for RPlillsi. 



VI. RELATION BETWEEN np AND A FREE 
ENERGY 

In this Section, the consequences of making the identi- 
fication np ^ (log Z) /N are explored. It has been shown 
some time ago that SP can be mapped onto the q ^ 1 
limit of the Potts model [43l li^. This mapping allows 
one to draw a parallel between an equilibrium thermody- 
namic transition (Potts) and percolation, which may be 
described as a purely geometric transition. One of the 
outcomes of this equivalence is the identification of the 
total number of clusters in percolation as the logarithm 
of the partition function for the resulting Potts model. 
The existence of this mapping has had profound impact 
on our understanding of scalar percolation, e.g. by allow- 
ing the exact calculation of percolation critical indices 
in two dimensions No equivalent mapping has been 
found for RP yet. However it has been suggested [20ll45| 
that a possible generalization for the free energy in RP 
for zero field is the total number of remaining degrees of 
freedom (zero modes or floppy modes). In the g — 1 case 
of RP each cluster has one remaining degree of freedom, 
thus the SP result is recovered. 

In Section Hvl it was shown that np (equivalently 7^) has 
to be continuous. This was used in Section^to locate the 
point 7c where the first-order transition happens. Thus 
np plays a role similar to that of the Gibbs free energy, 
which is also a continuous function of its intensive param- 
eters. In this section we add further evidence supporting 
the identification oinp as the logarithm of the "partition 
function" , or free energy, for the RP problem on random 
graphs. We demonstrate that the ansatz Zj^p oc e^"^ is 
consistent also in the presence of a field, and show how 
several thermodynamic quantities result from derivatives 
of this free energy. 

In Appendix lA II it is shown that the derivatives of np 
are 



dni 



57 

dnp 



dh 



i?- 1, 



(37a) 
(37b) 



As expected, the order parameter R results form deriving 
the logarithm of the partition function with respect to the 
field-like parameter, up to a constant shift. 
From Eqns. (|37|l we furthermore obtain 



dnf 



dnp 



1 



(i?2 - 1) + H{R-1) 



H 



dH 



7(i?-l). 



(38a) 



(38b) 
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A. Maxwell's Rule of Equal Areas 

The continuity requirement (|30|l employed in this work 
to locate the true transition point 7c is analogous to the 
condition that the free energy be continuous, as noted 
in previous work j20| |. The continuity of an appropriate 
thermodynamic potential is discussed in elementary sta- 
tistical mechanics textbooks in connection, for example, 
with Van der Waals' equation ^^|^, where it is used 
as a criterion to locate the transition. The continuity 
condition for our RP problem is written in general as 



dnp = 0, 



(39) 



where A and B are the two phases that coexist at the 
first-order transition (see Fig.EJ). Clearly H30(l is the par- 
ticular case of (|39|l for which the integration is done along 
a path of constant field. Choosing a path A ^ B on 
which 7 is constant instead, it is easy to put H39I) in the 
form of the famous Maxwell's "rule of equal areas" . Start 
from 







dh 



(40) 
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FIG. 7: This plot of the field H versus the order parameter R 
for Rigidity Percolation (g = 2 in this example) is the equiv- 
alent of a P-V diagram for fluids. Shown are the "isotherms" 
(solid lines) along which the coordination number 7 is con- 
stant. The thick dashed line delimits the "coexistence" region, 
and was obtained by solving H41|l numerically. The areas de- 
limited by the thin dashed lines above and below the (solid) 
horizontal coexistence line are equal. 



where the integral is done along an "isotherm" (a line of 
constant 7). After using (j37b|l and integrating by parts, 
this reads 



h{R,-f)dR, 



(41) 



whereji(i?, 7) = G-i(i?) 
Eq, 

rB 



l|lT|l is equivalent to the condition P{Vb — Va) = 
P{V,T)dV for a fluid, i.e. MaxweU's rule. Fig. □ 
serves to illustrate the fact that the areas above and be- 
low the horizontal coexistence line are equal, as H41|l re- 
quires. 



B. Energy, Entropy, and Work 

In the following we will assume that np (Eq. H29|l ) is 
the logarithm of the partition function, i.e. np — —7/ 
with / a free-energy density, and 7 = 1/T the inverse 
temperature. The energy density e is then 



dnp 



H 



\{l-E?) + H{l~R). (42) 



The entropy per site s turns out to be 

= np+^e^g[Pg{y)-{R-l)]. (43) 



dl 
dT 



H 



From the first law of thermodynamics, and using 
Eqns. H42|) and H43|) . we conclude that, in an infinitesimal 



transformation, the "work" done by the system against 
the environment is: 

dw = dq-de = -ds-de = -^dH = (R-l)dH (44) 
7 7 oH 

The analogous of the constant-pressure heat capacity for 
a fluid is in our case the constant-field heat capacity 



CH 



ds_ 



H 



2 9e 



H 



l-7Pg-i(2/)' 



(45) 



which is nonnegative by l|18|l . and diverges as I7 — 7*|~^ 
at the critical point. This of course should not be taken 
to mean that a — 1, since it is the constant-density heat 
capacity what defines a (For a discussion see Ref. _52J). 
The constant-field heat capacity cr is expected to diverge 
with the same exponent 7 as the susceptibility, unless the 
system has certain symmetries, which is the case for the 
Ising model, but not for RP. 

The analogous of the "constant-density" heat capacity is 



CR 



ds 



(46) 



which is zero and thus a = 0. 



C. Clausius-Clapeyron Equation 

Let A/ = -l/7AnF = -l/-i{np{B) - np{A)) be the 
free-energy difference between the two phases in the co- 
existence region. This is a function of the "temperature" 
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T = l/j and the field H 
texts ^3, we write 



/1/7. Following standard 



dH 
dT 



A/ 



dAf 




dT 


H 


dAf 




dH 


T 



where the derivative on the left-hand side is evaluated 
on any line of constant A/. Using 143|l one can write 
As = gi^Pg — AR). On the coexistence line, A/ = 
so Asj^^g^ — 7Ae|^^g^ and, using we have that 
^^Lex -f^'^R - 7/2A(i?2) = _Ai?(2^^±^). Thus 
the Clausius-Clapeyron equation for the RP problem can 
be written as 



dH 
'df 



VA 



(48) 



Since on approach to the critical point, y ^ y* — g ~ 1 
(Eq. (|13|l l. we conclude that the coexistence line reaches 
the critical point with slope (g — 1) in the H-T plane 
(See Fig. 13 . 



VII. DISCUSSION 

The Rigidity Percolation problem with g degrees of 
freedom per site has been considered on Erdos-Renyi 
graphs with average coordination 7. An external field 
h is introduced by connecting each site to a rigid back- 
ground, or "ghost site" , with a certain probability that 
depends on h as described in Section III CI The resulting 
Equation of State iQ for the density R of sites that are 
rigidly connected to the background, undergoes a first- 
order phase transition on a "coexistence line" "fdh)- This 
line ends at a critical point (Fig. |3J) with classical expo- 
nents a = 0, f3 — 1/2, 7=1 and 5 = 3. For comparison, 
the critical point of Scalar percolation (the case of 5 = 1) 
is located at h = and has different critical exponents: 
a = — 1,/3 = 1,7 = 1 and ^ = 2 Therefore in 
the MF approximation. Scalar and Rigidity percolation 
are in different universality classes. In two dimensions, a 
similar conclusion is reached by numerical means jl5ll45j |. 
were both transitions are continuous in zero field. 
It has been recently argued [s^ that for certain spin 
systems, the existence of a discontinuous transition in 
the MF approximation is enough to ensure that, in finite 
space dimensions d, there is a discontinuous transition if 
d is large enough. If a this result holds for RP, one would 
expect to see a discontinuous transition in zero field for d 
(and perhaps g) large enough. Up to now, extensive nu- 
merical simulations for RP have been performed only in 
two dimensions 0, |4^ ■ Except for patholo gica l 
cases where the transition happens at zero dilution |15|. 
it seems that the transition is continuous for all g in two 
dimensions. Large scale simulations are still needed to 
clarify this issue. 

On Cayley Trees with a boundary, g-rigidity percola- 
tion(RP) is closely related to {g + l)-bootstrap perco- 
lation(BP) (L9l i34]. Both percolate at the same density 



p of present bonds, although with different spanning clus- 
ter densities on trees. On random graphs, which have the 
local structure of a tree but no boundaries, fc-core tran- 
sition stays unchanged, but the RP transition is delayed 
to larger p values [l^- We have explicitely shown that, 
at the critical concentration 7c where the rigid cluster 
first appears, it has an exact balance of degrees of free- 
dom, i.e. it has exactly g bonds per site. This condition 
is satisfied both on the whole rigid cluster and on its 
{g + l)-rigid core. However globally the system has lesser 
bonds than needed to attain this balance, meaning that 
the rigid cluster is selectively made out of sites with more 
bonds than average. 

In previous work on Bcthc lattices in zero field [23| it 
was suggested that the number np of uncanceled de- 
grees of freedom is a free energy for the RP problem. 
For SP {g = 1) this result holds exactly, being one of 
the outcomes of the Fortuyn-Kasteleyn Random-Cluster 
model J^i, , which contains SP and the Potts model 
as particular cases. For RP, however, this identification 
only has the status of a plausible ansatz. This ansatz was 
used recently to predict some thermodynamic properties 
of chalcogenide glasses 56] , for which the RP transition 
has been shown to be relevant [llllT3 . IstIIs^ . 
Some of the reasons to believe that this identification 
might be correct in general are: a) the fact that np must 
be continuous at the discontinuous RP transition, and, 
b) the fact that for g = 1^ np is the number of connected 
clusters per site, so the FK result for SP is recovered ex- 
actly. In this work we have shown that, in the presence 
of a properly defined field /i, the ordered parameter R 
can be obtained as a derivative of the free energy with 
respect to h, thus adding further support to the belief 
that this identification is correct. 

Under the assumption that logZ = Nnp, the entropy 
per site s can be derived. The resulting expression 
was found to depend on the order parameter R alone. 
This is a property of other MF systems like for example 
the Ising ferromagnet. 

On the pedagogical side, we have attempted to situ- 
ate the discussion in terms of the parallel between the 
RP transition in a field and a condensation transition. 
Since only topological (connectivity) properties are im- 
portant for RP, it can be said that, in a sense, RP in 
a field is a sort of "geometric condensation transition" . 
Fig. [3 illustrates the similarities between both transi- 
tions. However a closer match is possible. Comparing 
{Tj) to lO one notices that, in a MF condensation transi- 
tion, A cx = 1/1;^ plays the role of the order parameter 
R in RP. From this point of view, the analog of a "vol- 
ume" in RP would be l/R^/"^. The analog of a P-V plot 
is shown in Fig. |S1 Notice however that, within this def- 
inition of a "volume" , P{V) does not satisfy Maxwell's 
rule of equal areas, V{P)dP = 0. 
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FIG. 8: Illustration of the similarities between RP and a 
condensation transition. In this plot, H — /1/7 is taken to be 
the analog of a pressure and the analog of a volume. 

The data are the same as shown in Fig. |7| 
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probability piP2, so that one can write 

{6} {8'}Ci3 

= J2Pp,iB)S{p,,B) (A3) 

{B} 

Notice that Pp^ {B') has a factor p2 for each present bond 
in configuration B\ and a factor l—p2 — 92 for each bond 
in B which is absent in B' . Factors associated with bonds 
not in B can be summed out to give one and therefore 
need not be considered. Thus if \B\ is the number of 
present bonds in configuration B and \B'\ the number in 

Now let P2 ^ 1- In this case we can expand S{p2,B) 
in powers of q2. The zeroth-order contribution comes 
from B' = B, the first-order contribution from the \B\ 
configurations B' which have exactly one bond b less than 
B, etc. 

S{p2,B) = pf + + 



beB 



2 \B\-2 
Q2P2 



J2 H'3~b-b') + ..., (A4) 



b.b'eB 



where we have written ;B — 6 to denote the result of delet- 
ing bond b from B. To first order in q2, 



S{p2,B) = TiB) + q2j2 iH13 -b)- , (A5) 

therefore 



beB 



APPENDIX A: DERIVATIVES OF THE NUMBER 
OF REDUNDANT BONDS 

Let J-{B) be an arbitrary function of the set of present 
edges B of a graph Q, and let P{B) be the probability to 
have edge set B. The average of J- is defined as 



{t) = J2p{b)HB), 

{B} 



(Al) 



where J2{b} ^ '-'^^^ configurations of present 
edges. 

If each edge is independently present with probability p 
and absent with probability 1 — p, then clearly 



Pp(S)-pl«l(l-p) 



A~\B\ 



(A2) 



where Bmax is the edge set of maximum possible cardi- 
nality. 

We now write p = P1P2 with < pi,P2 < 1- This can 
be realized by assuming that two types of bonds are in- 
dependently present with probabilities pi and p2 respec- 
tively, and that in order for an edge to be "active" , both 
types of bonds must be present on that edge. .F is now 
a function of active bonds only, which are present with 



{^)p.p. = {^)p. +<i2{Yl (-^(^ - ^) - •^(^)) ) ■ 

\beB / 

(A6) 

Taking derivatives with respect to p2 and letting P2 = 1 
one finally finds that, for arbitrary JF, 



dp 



(A7) 



\beB 



Eq. IIATII generalizes previous results of Coniglio in 
Ref. [531, and is the tool we use in this Section to de- 
rive some important relations. 

Consider now the case in which = Br{B), the total 
number of redundant constraints. In order to calculate 
the derivative of with respect to p, we notice that 
when removing a bond b from a configuration B of present 
bonds, the total number of redundant constraints will be 
reduced by one if and only if b is overconstrained. Oth- 
erwise if b is not overconstrained, Br remains unchanged 
by the removal of b. Eq. ljA7p then implies that 



d{B^_ d{Br)_ 
^ dp "^97 



(A8) 



where (-Sou) is the average number of overconstrained 
bonds in B. 
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A calculation of the field derivative of (Br) requires some 
additional considerations. Eq. HA7|I was derived under 
the assumption that each edge is independently present 
with some fixed probability, while our definition of ghost- 
field in Section Iff CI implies that each site can have any 
number n of ghost constraints with poissonian probabil- 
ity Pn{h) = e^^h^ /n\. We can represent this poissonian 
process in a way that enables us to use (jA7p as follows: 
we assume that each site has M » h slots connecting it 
to the background, and that each slot is occupied by one 
ghost edge with probability p/j = h/M . For large M, the 
probability to have n ghost edges on a site is Pn{h), so 
one has a realization of the poissonian distribution. This 
equivalent system has a total of NM ghost slots, each 
occupied by a ghost bond with a small probability ph , so 
we can now use Eq. (|A7|) . 

Eliminating a present ghost edge reduces the number of 
redundant constraints by one if both of the following con- 
ditions are satisfied: a) the remaining constraints that 
this site has are enough to rigidize (i.e. kj + kh > g) , and 
b) the remaining ghost-constraints that this site has are 
not enough to rigidize on their own (i.e. k^ < g). The 
joint probability for these two conditions to be satisfied 



g-i oo 

^p,(M-^p,(M J2 P'^i^)- 

j=0 j=0 k=0 

R{-/,h) - R{0,h) ^ R~ Rh, (A9) 

where we have defined Rh = Gg{h) as the density of 
rigid sites in the presence of the field h, for a graph with 
no edges (7 = 0). The total number of ghost slots is 
NM and thus the total number of present ghost bonds 
which upon removal produce a change in the number of 
redundant constraints is on average 



in agreement with H25() . and 



NM X — x(R- Rh) 
M ^ ' 

Eq. IjATfl . with ph — h/M, now implies that 

d{Br) 



dh 



N{R - Rh) 



Using {Br)-y — jrN/2 one then has 



dh 



2{R~Rh). 



(AlO) 



(All) 



(A12) 



(A14) 



in agreement with (|A12|I . Now since ^off — 5(1 ~ 
Gg+i{h)) — h[l — Gg{h)) one has that dgcs/dh = Gg{h) — 
1 — Rh — 1- Therefore the derivatives of np (Eq. (1^ ) 
are 



dnp 

dnp 
~dh 



i?2 



= R-1 



(A15) 



APPENDIX B: 



DERIVATION OF CRITICAL 
INDICES 



Eqns. H12af) and H12b|) imply that, for y y* , 



R = GM^R* + ^ + OiiAy)% (Bl) 

7' 



with Ay = y - y*, and y* ^ ^*R* + h*. The coefiicient 
of the 0{{Ay)^) term is nonzero. Furthermore Ay = 
Ah + 7* Ai?, + i?A7, and Eq. (|ljl|l can be rewritten as 

j*R ^ j*R* +x~x* +Ah + OiiAy)^) ^ 
RAj + Ah = 0{{Ay)^). (B2) 

On the critical isotherm (A7 — 0) we thus have that 
Ah - {j*AR + Sh)^ =^ AR^ (A/i)i/3, so (5 = 3. 
When h < h* the order parameter develops a 7-driven 
discontinuity Ai? oc (Ah)^ , which can be estimated in 
the following way. Assume R* is one of the three so- 
lutions of Eq. (|B1I) when Ah < 0, and determine the 
value of 7 for which this happens. Because Ay — and 
Ai? = 0, one finds that A7 = -l/R*Ah. Plugging this 
result back into Eq. ljBl|l and eliminating the solution 
R = R* one finds that the other two solutions behave as 
Ai?cx (A/i)i/2, or /3 = 1/2. 

The "susceptibility" x — ff U* diverges on approach to 
the critical point as x {Aj)^^ . 
gets 



Deriving Eq. ||7J) one 



dR 

dh 



1 " iPg-l 



(B3) 



1. Checking derivatives of 7r 

It is a trivial exercise to show that 7^ , as given by H28|) 
has the right derivatives, i.e. satisfies 



i?^ 



(A13) 



and recalling that 7* = {Pg-i{y*)) ^ we find that x cx 
(A7)-i, or 7 = 1. 

Notice that Rushbrooke's (a 2/3 7 = 2) and Griffith's 
{a + (3{S + 1)) — 2) relations are satisfied with a — 0, 
and this is consistent with the fact that the constant-i? 
specific heat (|46|) is zero. 
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